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1.  Introduction 


Asynchronism  in  the  arrival  of  power  at  the 
insulator-vacuum  interface  of  a  pulsed  power  machine 
leads  directly  to  a  loss  of  power  delivered  to  that 

interface'1,  in  addition  to  problems  associated  with 
possible  loss  of  insulation.  With  the  advent  of 
larger  pulsed  power  machines  composed  of  multiple 
modules  feeding  a  central  vacuum  region,  the  need  to 
reduce  the  asymmetry  in  the  power  density  delivered 
to  that  region  has  become  an  issue  of  some  importance. 
The  asymmetry  might  result  from  switching  in  the 
individual  modules  or,  as  we  discuss  here,  from  geo¬ 
metrical  effects.  In  this  paper  we  examine  how  such 

2 

geometrical  asymmetries  in  Double  EAGLE  might  be 
mitigated. 

2 

The  Double  EAGLE  pulsed  power  machine  at 
Physics  International  consists  of  two  transmission 
lines  which  feed  power  to  a  common  load.  Converting 
from  two  lines  to  a  central  axi symmetric  load  leads  to 
asymmetries  in  the  power  flux  density  around  the 
periphery  of  the  load,  as  is  shown  schematically  below. 


Figure  1.  Schematic  View  of  Double  EAGLE  Coupler 
Region. 

Consider  two  planar  strip  transmission  lines 
connected  to  a  circular  load.  The  power  pulse  will 
arrive  at  the  power-facing  points  of  the  arc  (left- 
right)  significantly  earlier  than  at  the  top-bottom 
points.  The  non-simultaneity  of  arrival  at  different 
points  on  the  load  (of  radius  R)  gives  rise  to  non¬ 
uniformity  in  the  voltage  circuit  and  hence  power. 

To  circumvent  this  problem  a  Transit  Time  Equali¬ 
zer  (TTE)  has  been  devised  to  increase  the  path  length 
for  signals  moving  near  the  center  of  the  strip-line. 

It  consists  of  a  hump  whose  dimensions  are  chosen  to 
add  a  delay  time  x  =  R/c  to  signals  propagating  over  it. 
An  artist's  impression  of  the  TTE  in  the  region  where 
coupling  of  the  linear  transmission  lines  to  the  cen¬ 
tral  load  takes  place  is  shown  in  Figure  2. 

2.  Theoretical  Model 

We  have  studied  computationally  the  propagation  of 
power  in  this  coupler.  Clearly  the  geometry  is  quite 
complex.  To  simplify  matters  we  have  examined  only  the 
propagation  of  signals  which  have  characteristic  wave¬ 
lengths,  A  =  ctQ  (where  t  is  the  characteristic  pulse 

time  and  c  the  signal  speed)  which  are  long  compared 
to  the  distance  A  between  the  opposing  conducting 


surfaces.  Under  this  simplification  the  electric 
field  E  is  constant  across  the  gap  or,  more  correctly, 
along  a  curve  that  intersects  normally  opposite  sur¬ 
faces.  We  need  consider  then  only  the  normal  compo¬ 
nent  of  the  E-field,  E,  and  may  ignore  all  variations 
in  the  normal  direction.  We  work  in  terms  of  the 


vector  potential.  A,  given  by 

*-■4  (1) 

-+  +  -> 

B  =  VxA  (2) 

V-A  =  0  (3) 

From  Maxwell's  equations 

a2 1  -* 

2-4-  =  -  cc  Vx  Vx  A  =  cT  A  (4) 

St^ 


Equation  4  is  strictly  valid  for  a  vacuum  line  or,  if 
we  replace  E  by  D  =  eE,  for  a  lirie  with  a  constant 

dielectric  e  such  that  c  =  c  £-  where  c  is  the 

o  o 

velocity  of  light  in  vacuum. 


Figure  2.  PFL  and  Coupler  Regions. 

To  solve  this  equation  we  use  the  usual  leap¬ 
frog  algorithm.  What  is  more  important  is  how  the 
discretization  of  this  equation  is  performed.  Briefly 
the  computational  mesh  shown  in  Figure  3  is  broken 
into  a  set  of  volume  elements.  Assuming  a  bilinear 
variation  of  the  field  quantities  E,A  over  these 

2 

finite  elements,  the  operator  V  can  be  constructed 
using  the  finite  element  method.  The  finite  element 
formulation  is  described  in  the  appendix. 
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Figure  3.  View  of  the  computational  mesh,  repre¬ 
senting  half  of  the  top  side  of  one  EAGLE 
module.  Using  mirror  symmetry,  this  is 
sufficient  to  represent  the  full  machine. 

3.  Calculations 

Figure  3  shows  a  cut-away  of  the  region  of 
Double  EAGLE  which  we  treat  -  namely,  the  region 
extending  from  the  pulse  forming  line  (PFL)  switch 
to  the  periphery  of  an  annular  resistive  load.  In 
the  work  presented  here  the  voltage  at  the  PFL  switch 
is  specified  (Figure  4).  For  the  purpose  of  this 
study  -  namely,  to  see  what  effect  the  TTE  has  on 
the  propagation  of  the  pulse,  -  this  is  adequate.  The 
dimensions  and  pulse-shape  correspond  to  specifica¬ 
tions  provided  by  Physics  International. 


The  inner  and  outer  surfaces  (cathode  and  anode) 
of  the  transmission  line  are  treated  as  perfectly 
conducting  and  the  water  medium  has  a  constant 
dielectric  in  which  the  signal  velocity  is 

c.  3  x  10^ e^2 m/s 
~  3.3  x  10^  m/s 

We  assume  that  A,e"  have  components  only  normal  to  the 
cathode-anode  surfaces,  and  that  mirror  symmetry 
exists  at  all  boundaries  other  than  the  PFL  input 

^specified  E)  and  the  load  arc  (matched  impedance), 
i hi s  is  exactly  correct  (assuming  both  modules  and 
all  switches  operate  in  perfect  synchrony)  except  for 
the  physical  edges  of  the  TTE  and  side-rails,  where 
it  is  a  good  approximation  and  corresponds  to  neglect¬ 
ing  perturbation  of  power  propagation  by  the  fringing 
fields. 


We  now  proceed  to  see  what  effect  the  TTE  has  on 
the  propagation  of  the  pulse  by  comparing  results  with 
those  obtained  in  a  geometry  where  the  TTE  is  removed. 
Figures  5 (a-h)  show  at  four  times  (t  =  99,  124,  149, 
and  223  ns)  the  voltage  levels  in  the  TTE  and  non-TTE 
geometries.  The  pulse  arrives  at  the  load  (Figs.  5c, d) 
with  a  front  which  more  closely  matches  the  circular 
load  in  the  TTE  geometry  than  in  the  non-TTE  case,  so 
that  vacuum  region  voltage  (and  power)  gradients  are 
initially  much  smaller  in  the  TTE  geometry.  There  is, 
however,  a  tendency  for  the  wave  to  form  itself  to  the 
load.  At  149  ns  (Figs.  5e,f)  the  flat  geometry  is 
slightly  superior  to  the  TTE.  Figures  6  and  7  sum¬ 
marize  the  load  voltage  asymmetries.  We  see  that 
during  the  initial  voltage  rise  the  TTE  geometry  does 
reduce  considerably  the  angular  variation  of  the  vol¬ 
tage  at  the  load.  Subsequent  "mode  sloshing"  results 
in  alternately  decreasing  and  increasing  in  the  voltage 
variations  in  both  geometries;  however,  the  amplitude 
tends  to  be  smaller  in  the  TTE  geometry. 


Figure  7.  Voltage  (flat  geometry)  as  a 
function  of  angle  at  the  load 
arc  for  various  times. 
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Figure  5a.  Pulse  in  the  TTE  geometry  at  99  ns. 

Peak  of  pulse  is  the  second-rightmost 
contour  band.  Note  delay  at  TTE. 


Figure  5c.  Pulse  in  TTE  geometry  at  124  ns.  Pulse 
is  relatively  uniformly  incident  on 
load  at  45  . 


Figure  5e.  Pulse  in  TTE  geometry  at  149  ns. 


Figure  5g.  Pulse  in  TTE  at  223  ns,  when  the  peak 
reaches  the  load. 


Figure  5b.  Same  as  5a  in  the  flat  geometry. 


Figure  5d.  Pulse  in  flat  geometry  at  124  ns. 

Substantial  voltage  gradients  appear 
on  load  area. 


Figure  5f.  Pulse  in  flat  geometry  at  149  ns. 

Flat  is  slightly  superior  to  TTE  at 
this  time. 


Figure  5h.  Same  as  5g  for  the  flat  geometry. 
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4.  Conclusions 


We  have  presented  an  electromagnetic  model  of  the 
coupler  region  in  Double  EAGLE.  The  model  is  suffi¬ 
ciently  general  that  it  is  applicable  to  arbitrary 
transmission  line  geometries  in  the  long  wavelength 
(A  »  A)  regime. 

Based  on  the  results  obtained  the  Transit  Time 
Equalizer  proposed  for  use  in  Double  EAGLE  does  reduce 
the  asymmetry  resulting  from  the  geometry  of  the  coup¬ 
ler  regions. 
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Since  A  has  only  a  normal  component,  A^ 


VxA 1 2  = 

VAn 
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and  the  contribution  to  from  a  given 
1  1 

Mi?  =  / /  VNi  ‘  VNj  D(y,v) 


3p 


element  is 


3r 

8y 


dydv  . 


To  complete  the  prescription  we  find  VN.  by  solving 
the  three  equations  1 
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Appendix 

Finite  Element  Prescription  for  Laplacian  Operator 
Our  treatment  is  based  on  the  identity 

d3r  A*  V2A  +  d3  r  (vxA)2  =  surface  terms 

Using  finite  element  discretization,  we  may  write  the 
above  as 

V,  A,  L. .  A  +  A,  M.  .  A.  =  surface  terms, 

I  ~  I  a*  I  J  ~U  ~  I  a*  I  J  ~  J 


n  •  VN^  =  0 


3r 


3r 

3v 


where  fi  =  -4^-x  •  Also,  the  contribution  of  an 


element  to  a  node  volume  is  given  by 
1  1 


(e)  = 


D(P,v) 


0  0 


sr 

3y 


8r 

"Sv" 


dy  dv 


where  Vi  is  the  volume  associated  with  each  node  of 
the  mesh.  We  then  identify  the  Laplacian  operator  as 


In  our  case,  we  have  a  quadrilateral  mesh  with 
the  interelectrode  spacing,  D,  and  one  component  of  A 
defined  at  each  node.  The  matrix  M^.  is  formed  as 

the  direct  sum  of  contributions  from  each  quadri¬ 
lateral  element  as  follows: 


Define  a  transformation  T  between  the  quadri¬ 
lateral  and  the  unit  square  as 

r  =  T(P,v)  =  (l  -  y) ( 1  -v)rx  +  y(i  .  v)r2  +  yvr3 


+  (1  -  y)v  r4 
4 

=  2  My,v)r. 
i  =  1 


where  ^...r^  are  the  actual  coordinates  of  the 
corners  of  the  quadrilateral.  We  interpolate 
4 

A(y,v)  N.(y,v)Ai 

i  =  1 

4 

D(y»v)  N.  (y,v)D. 

i  =  1 
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